Simultaneous Spin-Charge Relaxation in Double Quantum Dots 
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We investigate phonon-induced spin and charge relaxation mediated by spin-orbit and hyper- 
fine interactions for a single electron confined within a double quantum dot. A simple toy model 
incorporating both direct decay to the ground state of the double dot and indirect decay via an 
intermediate excited state yields an electron spin relaxation rate that varies non-monotonically with 
the detuning between the dots. We confirm this model with experiments performed on a GaAs 
double dot, demonstrating that the relaxation rate exhibits the expected detuning dependence and 
can be electrically tuned over several orders of magnitude. Our analysis suggests that spin-orbit 
mediated relaxation via phonons serves as the dominant mechanism through which the double-dot 
electron spin-flip rate varies with detuning. 



Controlling individual spins is central to spin-based 
quantum information processing [1-3] and also enables 
precision metrology [ , ]. While rapid control can be 
achieved by coupling the spins of electrons in semiconduc- 
tor quantum dots [ , 2] to electric fields via the electronic 
charge state [ , 6-14], spin-charge coupling also leads 
to relaxation of the spins through fluctuations in their 
electrostatic environment. Phonons serve as an inherent 
source of fluctuating electric fields in quantum dots 
and give rise to both charge and spin relaxation through 
the electron-phonon interaction. In GaAs quantum dots, 
the direct coupling of spin to the strain field produced by 
phonons is expected to be inefficient [' ^ , ^^]. The dom- 
inant mechanisms of phonon-induced spin relaxation are 
therefore indirect and involve spin-charge coupling due to 
primarily spin-orbit [^ H-9i] and hyperfine ['^O-O't] inter- 
actions. Confining an electron within a double quantum 
dot (DQD) provides a high degree of control over the 
charge state ['^^-'^^], so that relaxation rates can be var- 
ied over multiple orders of magnitude by adjusting the 
energy level detuning between the dots [25, 32-35]. 

Here, we investigate the interplay of spin and charge 
relaxation via phonons for a single electron confined to 
a DQD in the presence of spin-orbit and hyperfine in- 
teractions. We present a simple model together with 
measurements of the electron spin relaxation rate in a 
GaAs DQD, both of which yield a non-monotonic depen- 
dence on the detuning between the dots (see Fig. 3). 
The experiments provide confirmation of the model and 
demonstrate the existence of spin "hot spot" features [ 8- 
20, 3(3, ^ ] at nonzero values of detuning, where relax- 
ation is enhanced by several orders of magnitude. The 
opposite behavior is observed at zero detuning, where the 
spin-flip rate exhibits a local minimum. Theoretically, 
spin hot spots are predicted to occur due to the com- 
plete mixing of spin and orbital states at avoided energy 
crossings associated with spin-orbit coupling [ o, ly, o ]. 
From a practical standpoint, adjusting the detuning to 
these points represents a potential method for rapid all- 
electrical spin initialization. 



We describe a single electron confined within a DQD 
(Fig. 1) using a toy model that includes only the lowest- 
energy orbital level of each dot. This two- level approxi- 
mation [ ] enables the charge degrees of freedom to be 
represented by the Pauli matrices , a^y , in the basis 
, where \L) {\R)) denotes an electron in the left 
(right) quantum dot and = \L) {L\ — \R) {R\. We can 
then express the Hamiltonian of the system as 

Hd = i^O + ^so + ^nuc, (1) 

^0 = -^o-z - tax - ^zSz, (2) 

^so = Kso • S cr^, (3) 
^nuc = Knuc ' S (J^. (4) 

The first two terms in Hq [Eq. (2)] specify the purely or- 
bital part i^orb = (e/2) (j^ — tax of the electronic Hamil- 
tonian in terms of the energy level detuning e and the 



Empty dot 




wait for time T 




Figure 1: Electrochemical potential diagrams for a dou- 
ble quantum dot (DQD), illustrating the measurement cycle 
used to obtain the experimental spin relaxation rate (see main 
text). Varying the detuning e between the left (L) and right 
(R) dots while keeping the tunnel coupling t fixed (stage 3) 
tunes the relative energies of the charge states. Tunneling of 
the electron between the dots is accompanied by spin rotation. 



tunnel coupling t between the two dots (Fig. 1). Diago- 
nalization of i^orb yields the eigenstates 

1+) ^ cos^|L)-sin^|i?), 

I-) ^ sin^|i:)+cos^|i?), (5) 

where (j) varies with e and t according to tan(/) = 2t/e. 
The corresponding eigenvalues are separated in energy 
by a gap A = - E_ = + 4^^ [see also Fig. 2(a)]. 

Spin dependence is introduced into the Hamiltonian 
via the last term in i^o, together with H^o and i^nuc [see 
Eqs. (2)-(4)], where the vector of electron spin oper- 
ators is denoted by S = {Sx^ Sy^ Sz)- The term H^o 
describes spin-orbit coupling which is linear in the elec- 
tron momentum p. The general form given in Eq. (3) 
takes into account both the Rashba [38, 39] and the lin- 
ear Dresselhaus [ ] forms of spin-orbit interaction, with 
strengths and orientations that are specified by the vec- 
tor Kso = (r, Note that H^o acts as ay within 
the orbital subspace, which follows from parity selection 
rules for the matrix elements of p in the basis , 
Thus, Hso describes tunneling between the dots accom- 
panied by a spin flip (Fig. 1) [ ]. 

The remaining spin-dependent terms in represent 
forms of the Zeeman interaction that are distinguished by 
their action within the orbital subspace. The final term 
in Hq represents the coupling of the electron spin to a 
magnetic field of strength B = Qz/ \9\ I-(^b that is uniform 
over the two dots, where g is the electron g-factor and 
/is is the Bohr magneton. The vector Knuc = {u^v^w) 
in i^nuc specifies the strength and orientation of a mag- 
netic field gradient across the two dots, i^nuc acts as 
within the orbital subspace. For GaAs quantum dots, 
^nuc can be used to model the hyperfine interaction be- 
tween the electron spin and the ensemble of nuclear spins 
with which the electron wavefunction overlaps. The as- 
sociated intrinsic magnetic field gradient is assumed to 
originate from an effective nuclear field Bnuc with a ran- 
dom, spatially-varying orientation described by a Gaus- 
sian distribution and magnitude 5nuc given by the root- 
mean-square (rms) value [23-25, 27, 4 ]. 

Figure 1 illustrates the scheme used for the measure- 
ment of the spin relaxation rate. The experimental setup 
is described in [42]. In the first step of the measurement 
cycle, a single electron spin is initialized by emptying the 
DQD and then letting a single electron tunnel into the 
left dot far from the degeneracy of \L) and \R). The elec- 
tron spin is randomly up or down. Next, a voltage pulse 
adjusts the electrochemical potential of the right dot to 
tune the detuning closer to the degeneracy to a value e 
for a wait time r. After the wait time, the electrochemi- 
cal potential is pulsed back and the spin of the electron is 
read out using energy- selective spin-to-charge conversion 
[43]. This cycle is repeated for a given e and r to obtain 
an average spin-down probability at the end of the cycle. 



For each series of measurements as a function of r at a 
fixed e, the amount of detected spin-down is fitted with 
an exponential decay, from which the spin-relaxation rate 
at each e is obtained as shown in Fig. 3. 

The variation of the measured spin relaxation rate with 
detuning can be understood in terms of the spectrum for 
the one-electron double dot. Figure 2(a) shows an exam- 
ple spectrum for Hd [Eq. (1)] as a function of detuning. 
In Fig. 2 and throughout the present work, we consider 
the limit t <^ Qz which corresponds to the measurements 
described above (see [4^^]). The notation t, i used to label 
the states in the figure refers to the components of spin 
along the quantization axis defined by the external mag- 
netic field. In accordance with the experiment [ ^'^, 44]^ 
we choose this field to be in the plane of the quantum dots 
and parallel to the double-dot axis. The in-plane crys- 
tal lattice orientation characterizing the spin-orbit inter- 
action [Eq. (3)] is parametrized by an angle relative 
to this axis. Of particular consequence for the present 
work is the fact that Hso gives rise to avoided crossings 
in the spectrum at e ~ i^^, where maximal coupling 
of the states |+,t) and | — ,i) occurs and leads to the 
complete mixing of orbital and spin degrees of freedom. 
These finite values of e correspond to spin "hot spots" 
[ )-20, 36, 37] and are associated with enhanced spin 
relaxation rates, as is shown below. 

Including coupling to phonons in the description of the 
single-electron double-dot system leads to the Hamilto- 
nian H = Hd + i^ep , where 




X (a.,k + <_k) e^^" (6) 

is the electron-phonon interaction [ ], expressed in 
terms of the mass density po, the volume Vq, the phonon 
speeds c^^, the deformation potential and the piezo- 
electric constant S. The operator aj, ^ (a^^^k) creates 
(annihilates) a phonon with wavevector k and polar- 
ization u [the sum over u is taken over one longitudi- 
nal (/) mode and two transverse {t) modes], and S^^i 
is the Kronecker delta function. Fermi's golden rule 
for the rate F of phonon-induced relaxation of the elec- 
tron from state \i) to state |/) of the double dot gives 
F ~ |(/| i^ep P {^d)i where p (A^) is the phonon den- 
sity of states evaluated at the gap A^^ between levels i and 
/ that determines the energy of the emitted phonon. 

We first consider a qualitative model for F, where we 
estimate the transition matrix element (/| e^^'^ \i) (see 
[ ]) by writing e^^'^ ~ 1+zk-r and determining the corre- 
sponding matrix element of the dipole operator d = — er 
(here, e denotes the magnitude of the electron charge). 
To evaluate dipole matrix elements, we define Gaussian 
wavefunctions (r) = (r| L (R)) which are shifted 

along the dot axis by ±a/2 for the left-localized and right- 
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localized orbital states. While ijjL and ijjR are not orthog- 
onal, their overlap is small. We neglect corrections due 
to this overlap in our calculations. Using these wavefunc- 
tions, we find d = Dz with D = (ea/2) cr^. The qualita- 
tive form of the relaxation rate can then be approximated 
by r ~ \df F (Ad) , where d denotes the first-order term 
of {f\ D\i) and F {A^) represents the dependence of the 
rate on the gap energy A^^ (see l^"^] for more details). 

To identify the states of the double dot between which 
phonon-induced relaxation occurs, we treat V = H^o + 
^nuc [Eqs. (3) and (4)] as a perturbation with respect to 
Ho [Eq. (2)] and use nondegenerate perturbation theory 
(which is valid away from e ~ i^^) to find the first-order 
corrections to the energies and eigenstates of Hq. We de- 
note the corrected states by || — , t), | — |+, t)^ 
and consider relaxation of the electron spin from the 
excited state | — ,4) to the ground state | — ,t) of the 
DQD [see Fig. 2(a)], which can occur directly as well 
as indirectly via the excited state |+,t)- Away from the 
avoided crossing points, we note that |+,t) ^ l+^t) and 
| — ^ | — ,t)- The state |+,t) therefore relaxes rapidly 
to | — ,1)7 as effectively only orbital decay is involved and 
no spin fiip is required in this second step [46]. In the 
following, we assume that this charge relaxation is in- 
stantaneous compared to the spin relaxation and use the 
dipole matrix element for | — ,i) l+^t) to describe the 
full indirect transition. 

We approximate the relaxation rates and Fg [Fig. 
2(a)] in the presence of both H^o and i^nuc by calculat- 
ing the first-order terms and d^ of the dipole matrix 
elements (— — ,^) and (+, — , ^), respectively. 
These terms are functions of the spin-fiipping compo- 
nents r, s, li, and v in Eqs. (3) and (4). Averaging over 
the nuclear field distribution [23-25, 27] 
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where 



^nuc = I^I^S^nuc = y(|Knuc| )/3, givCS (u) 

{v) = and (u^) = (v^) = hl^^. We thus find 



{\de\') = 



ea ( 2t 



2 \AJ {A-n^){A^n^) 
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2 I 2 

r + s 



2t 



X, (8) 
(9) 



(26; 



nuc) 



These expressions are plotted in Fig. 2(b). Note that 
both Eqs. (8) and (9) are undefined at the avoided cross- 
ing points in Fig. 2(a), where A = Qz- Thus, the curves 
shown in Fig. 2(b) are valid only away from these points 
(i.e., where nondegenerate perturbation theory is a rea- 
sonable approximation) . Both (|<i6p) and (|<ieP) are only 
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Figure 2: (a) Spectrum of Hd [Eq. (1)] as a function of 
detuning e for the case t ^ r^^, in the presence of spin-orbit 
coupling (Kso 7^ 0, Knuc = 0). Correspondence with the 
eigenstates of Hq [Eqs. (2) and (5)] is indicated for relevant 
regions of the spectrum. Avoided crossings due to spin-orbit 
coupling occur aX e ^ ±Qz- The spectrum shown corresponds 
to t = 10 /ieV [42], B = 6.5 T [42, 44], the dot size a = 15 nm 
and interdot separation a = 110 nm, the GaAs effective mass 
m* = 0.067me (where rUe is the free-electron mass) and g- 
f actor g = 0.36, the Rashba and linear Dresselhaus spin-orbit 
strengths ao = 3.3 x 10~^^ eV-m and /3o = 4.5 x 10~^^ eV-m, 
respectively, and = 7v/8. The values of ao, /^o, and are used 
to determine Kso = (r, s,^). (b) Dipole-dependent factors 
and {\de\'^) [Eqs. (8) and (9)] used to qualitatively 
model the relaxation rates Tb and Fe in (a), as a function of 
detuning for Bnuc — (solid lines) and Bnuc — 3 mT (dotted 
(7) lines), with ao = 3.3 x 10"^^ eV-m and /3o = 4.5 x 10"^^ eV-m. 
All other parameters are identical to those used in (a). Units 
for the dipole are given in terms of the Bohr radius ao- The 
dipole model is not valid in the shaded region, (c) Relaxation 
rates Fa, F?,, and Fe as a function of detuning. Fhe rates are 
calculated using pQ — 5.3 x 10^ kg/m^, ci — 5.3 x 10^ m/s, 
ct = 2.5 X 10^ m/s, A = 7.0 eV, and S = 1.4 x 10^ eV/m 
[ ], together with the parameter values used in (b). Lines 
are guides to the eye. 



slightly modified by the coupling of the electron spin to 
an effective nuclear field of rms strength ^nuc = 3 mT 
['^^], as expected from Eqs. (8) and (9) in which the 
nuclear field term is scaled with respect to the spin-orbit 
terms by a factor (21/0.^)^ <^1[ ]. Saturation of (|<i6p) 
occurs at zero detuning for both the 5nuc = and the 
^nuc = 3 mT cases. On the other hand, (|<ieP) vanishes 
at e = regardless of the value of B^^^- As the experi- 
mental relaxation rate contains a local minimum at zero 
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Figure 3: Experimental detuning-dependent single-spin re- 
laxation rate (Fexpt) and comparison with the toy model de- 
scribed in the present work (Fth) for both zero and finite 
temperature. Error bars depict 90% confidence intervals for 
the data. The parameter values used to calculate Fth are the 
same as those used in (b) and (c) of Fig. 2. 



detuning (see Fig. 3), the present analysis suggests that 
the direct transition | — ,i) \ — alone does not ac- 
count for the observed spin relaxation and that indirect 
relaxation via the excited state |+,t) potentially plays 
a significant role in the spin- flip rate. The relative con- 
tributions of the direct and indirect transitions to the 
overall rate are explored in [^'•]. 

To compare our theoretical predictions more directly 
with the experimental results, we carry out the full cal- 
culation of the relaxation rates for both direct and indi- 
rect transitions to the ground state by applying Fermi's 
golden rule to relaxation induced by i^ep [Eq. (6)]. De- 
tails are given in the Supplemental Material [ ] . We set 
Knuc = for simplicity, as the preceding analysis based 
on dipole matrix elements suggests that the hyperfine 
term i^nuc represents only a small correction to the de- 
cay rate [see Eqs. (8), (9) and Fig. 2(b)]. Application 
of a Schrieffer- Wolff transformation [ ] enables diago- 
nalization of the full double-dot Hamiltonian for all 
6, including the avoided crossing points e ^ i^^, and 
the eigenstates of Hd are used to calculate the relaxation 
rates via Eqs. (S1)-(S3) [ ]. 

Relevant portions of the curves for the decay rates r2i, 
F31, and F32 (where we number the levels according to 
their energy eigenvalues and use F^ f to denote the rate of 
relaxation from level i to level /) are plotted together in 
Fig. 2(c). The rate F^ is associated with mainly charge 
relaxation and is given by F21 (Fsi) for |e| < ftz (|e| ^ 
r^^), while F5 is associated with mainly spin relaxation 
and is given by F31 (F21) for |e| < (|e| > Q,). The 
rate Fe corresponds to a combination of spin and charge 
relaxation and is given by F32 for all |e|. Note that F^ 
Fg, which is consistent with our prior assumption that 
the effective rate for indirect relaxation to the ground 



state is determined by Fg. 

For |e| ^^zi indirect spin relaxation occurs by a tran- 
sition to the lower-energy intermediate state via phonon 
emission [Fig. 2(a)]. On the other hand, the indirect 
transition for |e| > Qz requires phonon absorption in 
order to excite the electron to the higher-energy inter- 
mediate state. Using the Einstein coefficients and the 
Bose-Einstein distribution (n) = 1/ [exp (Ac^/Zc^T) — 1] 
(where /c^ is the Boltzmann constant and T is the tem- 
perature) to express the rates of spontaneous emission, 
stimulated emission, and absorption associated with the 
lowest three double-dot levels in Fig. 2(a) in terms of 
Fq, F^,, and Fe [32], we take the full theoretical detuning- 
dependent spin relaxation rate Fth to be given by F5 + Fg 
for |e| < Qz and by F5 + Fg (n) / ((n) + 1) for |e| > Q^- 
This rate is plotted together with the zero-temperature 
limit of the model and the measured rate Fexpt in Fig. 3 
for T = 250 mK [42, 44]. At the avoided crossings asso- 
ciated with spin-orbit coupling (e ^ ±^2)7 we find peaks 
in Fth that resemble the spin hot spot peaks observed ex- 
perimentally. The zero-detuning minimum found in the 
measurements appears in both the zero-and the finite- 
temperature models. In addition, close qualitative agree- 
ment between the finite-temperature model and experi- 
ment is observed for a wide range of detuning values. 
While limitations of our theoretical description arise from 
the two- level approximation we use for the orbital states, 
we have nevertheless shown that several characteristic 
features present in the measured detuning dependence of 
the double-dot spin relaxation rate can be understood 
within this simple model. 

The results of the present work therefore suggest that, 
in accordance with the case of single lateral GaAs quan- 
tum dots [ ], the observed variation of the spin relax- 
ation rate with detuning in double dots is dominated by 
spin-orbit mediated electron-phonon coupling. The spin- 
orbit interaction may thus provide the key to rapid all- 
electrical initialization of single spins. 

Note added. During the preparation of this 
manuscript, we became aware of a recent experimental 
observation of a spin hot spot in a Si quantum dot [ *8]. 
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I. SPIN RELAXATION RATE MEASUREMENT 

The measurement illustrated in Fig. 1 of the main text 
is performed on a DQD formed by using Ti/Au surface 
gates to locally deplete a two-dimensional electron gas 90 
nm below the surface of a GaAs/(Al,Ga)As heterostruc- 
ture. Quantum point contacts (QPCs) on both sides of 
the DQD are used to detect the charge on each dot. We 
tune the DQD to contain a single electron (Fig. SI) and 
operate close to the degeneracy of \L) and \R). To set 
the electrochemical potentials in the left and right dot in- 
dependently, combinations of voltages on two gates, one 
closer to the left dot and the other closer to the right 
dot, are tuned which compensate for capacitive cross- 
coupling. The tunnel coupling is tuned to be approxi- 
mately 8 jneV. Both the tunnel coupling and the con- 
version from applied gate voltage to change in detuning 
are determined from microwave spectroscopy [ ]. An in- 
plane magnetic field ^ext = 6.5 T is applied to split the 
spin-up and spin-down energy of the electron by the Zee- 
man energy {Ez ~ 130 jaeV). The electron temperature 
is typically 250 mK [_]. 
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Figure SI: Schematic of the double quantum dot charge sta- 
bility diagram associated with the steps in the measurement 
scheme of Fig. 1 in the main text, (n, m) represents a state 
with n electrons in the left dot and m electrons in the right 
dot. The charge configurations (1, 0) and (0, 1) correspond to 
the states \L) and \R), respectively, introduced in the main 
text. AVl and AVr represent the effective voltages applied 
to the left and right gates. 



II. COMPARISON OF DIRECT AND INDIRECT 
SPIN RELAXATION RATES 



The relative contributions of the direct and indirect 
transitions described in the main text to the overall re- 
laxation rate are characterized by the ratio Fe/F^. An 
estimate of this ratio can be obtained by combining 
the dipole-dependent factors (|<i6p) and (|<ieP) given in 
Eqs. (8) and (9) of the main text with the correspond- 
ing energy-dependent factors F (A5) and F (Ae) , where 
A5 ^ and Ag ^ Qz — ^ are the associated energy 
gaps. To determine an approximate form for F (A^^), we 
consider the full expression obtained from Fermi's golden 
rule for the rate of decay induced by i^ep [Eq. (6) of the 
main text]. Defining = e*^ '", we obtain 
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S^) hng{k)S (Sph - Ad) , 



(SI) 



where we have converted the sum over k to an integral, 
and 



Lnsik) = J dn\{f\Mk\i)\' 



(S2) 



is a momentum-space angular integral that is determined 
by the transition matrix element (/| In Hq^ and 

Eq. (SI), a linear dispersion = hcj^k for the phonon 
energies is assumed. The only factor in i^ep which de- 
pends on electronic degrees of freedom (specifically, on 
the position operator r for the electron) is M^, which 
can be expressed within the same two-level approxima- 
tion used for Hd in Eqs. (l)-(4) of the main text. The 
form of Mk in the basis {\L) ,\R)} is given by 
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akz 



- zsm 1 -- I cr^ + e '''/^^V^ 



, (S3) 



where the term involving cFx describes phonon-assisted 
tunneling between the left and right dots. 
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Figure S2: Comparison of the ratio Fe/Fb of the relaxation 
rates plotted in Fig. 2(c) of the main text (points) with the 
approximation in Eq. (S4) obtained using dipole matrix ele- 
ments (line). 

Retaining only the piezoelectric term (the dominant 
term for the parameters we choose [ ]) in Eq. (SI), set- 
ting = = hcj^k^ and noting that /ang(^) ^ l/k'^ 
in the limit of large [ ] , which is the limit appropriate 
for the size of A^^ at small detuning |e| [see Fig. 2(a) of 
the main text], we find F {Ad) ~ ^d^- The ratio of the 
decay rates [averaged over the nuclear field distribution 



in Eq. (7) of the main text] can then be estimated as 




As e -> 0, A ^ 2t <C SO that (re)/(r5) approaches 
(e/2t)^ . The estimated ratio in Eq. (S4) agrees well with 
the ratio re/r5 obtained from Fig. 2(c) of the main text 
for |e| < (Fig. S2). 
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